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Abstract. Inclusive cross-sections for gluon jet production are studied numerically in the 
perturbative QCD pomeron model for pA and central AA collisions at high energies. Two forms 
for the inclusive cross-sections, with and without emission from the triple pomeron vertex, are 
compared. The difference was found to reduce to a numerical factor ~ 0.7 -j- 0.8 for momenta 
below the saturation momentum Q s . Above Q s no difference was found at all. For pA collisions 
the gluon spectrum was found to be ~ A 7 at momenta k below Q s and ~ A ' 9 above it. For 
central AA collisions it was found to be ~ A at momenta k below Q s and ~ A 1, above it. At 
large k the spectrum goes like 1/fc 2 ' 7 ^ 3,3 flattening with energy. The multiplicities turned out to 
be proportional to A 0,7 for pA collisions and A for central AA collisions with a good precision. 
In the latter case they are becoming more peaked at the center with the growth of energy. Their 
absolute values are high and grow rapidly with energy in accordance with the high value of the 
BFKL intercept. 



1 Introduction 

New experimental data on heavy-ion collisions at RHIC and future such data to be obtained 
at LHC attracted much attention to the spectra of produced secondaries and in particular 
their dependence on the atomic number of colliding nuclei. One would like to have relevant 
predictions based on the fundamental theory and not purely phenomenological. At present the 
only candidate for this is the hard pomeron model derived from perturbative QCD. Originally 
constructed for the description of high-energy low -x hadronic scattering (the BFKL model [1] ) 
it has subsequently been generalized to hadronic or deep inelastic scattering on nuclei [2, 3] 
and nucleus-nucleus scattering [4]. The model suffers from a serious drawback related to the 
use of a fixed and not running strong coupling constant. Curing it does not look too promising, 
since due to absence of ordering of momenta in the model, it also means solving the confinement 
problem. However in spite of this defect the model seems to describe high-energy phenomena 
in a qualitatively reasonable manner. Also attempts to include the running of the coupling in 
some effective way have shown that the effect of the running is not at all overwhelming, although 
introduces some quantitative changes into the predictions. So, also for lack of something better, 
the perturbative QCD pomeron model appears to give a reasonable basis for the discussion of 
particle production in very high-energy heavy-ion collisions. Of course due to the perturbative 
character of the model it can only give predictions for production of jets, leaving jet-to-hadrons 
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conversion to non-perturbative fragmentation mechanism. Also it has to be stressed that the 
model is, by construction, oriented towards very high (asymptotic) energies or, equivalently, very 
low values of x. So its application to present-day energies does not seem to be fully justified, 
at least in the lowest order of perturbation expansion to be used in the present calculations. 
This has to be kept in mind when comparing the prediction of the model with the existing 
experimental data. 

From the start it has to be noted that in the model scattering amplitudes can be found 
comparatively easily only for hadron-nucleus collisions. They are given by a set of all pomeron 
fan diagrams, which are summed by the non- linear evolution equation of [2, 3]. Nucleus-nucleus 
collision amplitudes are described in the model by complicated equations, whose solution is quite 
difficult to obtain even numerically (see [5] for partial results). Happily, as was shown in [6], 
due to Abramovsky-Gribov-Kancheli (AGK) cancellations [7], to find the single inclusive distri- 
butions one does not have to solve these equations, but only to sum the appropriate sets of fan 
diagrams, that is to solve the much simpler hadron-nucleus problem. Still the latter problem 
involves a numerical study of considerable complexity. So up to now there has been no consistent 
calculation of the jet spectra for realistic nuclei at very high energies, although some preliminary 
attempts has been done in [8, 9, 10, 11]. In all cases however the authors relied on very drastic 
simplifications, from the start choosing for the nuclear structure and/or for the gluon distribu- 
tions in the colliding nuclei some primitive explicit forms in accordance with their own taste 
and prejudice. In fact these forms appear to be rather far from realistic ones, which correspond 
to actual participants and follow from the calculations. This gave us motivation to calculate 
numerically the jet spectra in hadron-nucleus and nucleus-nucleus collisions as predicted by the 
hard pomeron model in a consistent manner. 

Another goal of the present calculations has been to compare the results obtained on the ba- 
sis of the expression for the inclusive cross-section which follows from the AGK rules applied to 
the diagrams with QCD pomerons interacting via the three-pomeron coupling [6] with a some- 
what different expression obtained from the colour dipole picture [12]. Our calculations show 
that these two formally different expressions lead to completely identical results at momenta 
of the order or higher than the value of the so-called saturation momentum Q s . At momenta 
substantially lower than Q s the colour dipole cross-sections differ from the ones from the AGK 
rules by a universal constant factor ~ 0.7 -j- 0.8. 

In both cases the yl-dependence of the spectra at momenta below Q s is found to be suppressed 
as compared to the naive probabilistic expectations, which predict that it should follow the 
number of binary collisions, oc A for pA scatering and oc A 4 / 3 for AA scattering at fixed impact 
parameter. Insted we have found a behaviour damped by roughly factor A 1 / 3 , that is oc A ' 7 for 
pA collisions and oc A for central AA collisions. Since Q s grows with energy very fast, the region 
where the spectra behave in this manner extends with energy to include all momenta of interest. 
At momenta greater than Q s the spectra grow faster but still much slowlier than the number of 
binary collisions (a numerical fit gives something like oc A 1 ' 1 for central A A collisions). 

Note that in the last years a few more phenomenologically oriented studies of particle produc- 
tion in nucleus- nucleus production at very high energies have been presented, in the framework 
of the color-condensate model [13] solved in the classical approximation on the lattice [10] 
and in the saturation model [11]. In both approaches quantum evolution of the nuclear gluon 
density was neglected and the saturation momentum was introduced as a parameter fitted to 
the experimental data at RHIC. Although some of their predictions (proportionality of the mul- 
tiplicity to A modulo logarithms) agree with our calculations with full quantum evolutions, the 
quantitative results are rather different. We postpone a more detailed discussion of this point 
until our Conclusions. 

The results of this study continue and extend the ones published in [14]. As compared to 
that publication we explore the y-dependence of the cross-sections and, apart from AA central 
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collisions, study pA collisons as well. We also concentrate on energies accessible in the near future 
(at LHC) to be closer to the experimental situation. Still we have to stress that our results do not 
pretend to describe the existing data obtained at RHIC. As mentioned, the perturbative QCD 
pomeron model studies phenomena at very high energies, at which the coupling constant becomes 
sufficiently small. The characteristic energetic variable for the model is the scaled rapidity 
Y = (N c a s /ir)Y. One observes that the solution of the non-linear evolution equation acquires 
its standard scaling form, independent of the choice of the initial conditions, at Y > 1.5 -j- 2. 
With a s ~ 0.2 this sets the lower limit for the rapidity at the center, at which our results can 
be applied, to be 7.54- 10, which implies the overall rapidity for the collision in the region 15-i- 
20. This is much higher than the rapidities at RHIC but well within the possibilities of LHC. 
So we may hope that our predictions can be tested at LHC but have little chance of success for 
the data from RHIC. In relation to the latter the model can possibly indicate some trends of the 
observable cross-sections with the growth of energy, much in the same manner as it predicted the 
growth of the total cross-sections with energy at the time of its appearance (although the rate of 
this growth was grossly overestimated). The found proportionality of the spectra to the number 
of participants and not to the number of binary collisions well agrees with the experimental 
findings and is one of these trends. 



2 Basic equations 
2.1 Total cross-sections 

As mentioned in the Introduction, the scattering amplitude can be relatively easily found for 
hA collisions but not for AB collisions. Corespondingly here we present formulas which serve 
to calculate the total cross-sections for hadron (proton)-nucleus collisions. This cross-section is 
given by an integral over the impact parameter b and pomeron transverse dimension r as 

a A {y) = 2 J d 2 bd 2 rp(r)$(y, r, b). (1) 

Here y is the overall rapidity, p(r) is the colour dipole density of the projectile (proton) and 
<J>(y, r, b) is the colour dipole-nucleus cross-section at fixed y, r and b, given by the sum of all fan 
diagrams stretched between the projectile and nucleus. 

Function <f>(y, r, b) = $>(y, r, b)/(2irr 2 ), in the momentum space, satisfies the well-known non- 
linear equation [2, 3] 

M|^ = _^(^, 6W W), (2) 

where y = ay, a = a s N c /ir, a s and N c are the strong coupling constant and the number of 
colours, respectively, and H is the BFKL Hamiltonian. Eq. ( 2 ) has to be solved with an initial 
condition at y = determined by the colour dipole distribution in the proton smeared by the 
profile function of the nucleus. In our calculations we take this distribution in accordance with 
the Golec-Biernat distribution [15], duly generalized for the nucleus: 

<P(0,Q,b) = ^AT A (b)a f(q), (3) 

where 

1 q 2 

f(q) = Ei(-x), x = ± 7 , (4) 

W 2 v h 0.218 GeV 2 

ao = 20.8 mb and T A {b) is the standard nuclear profile function, which we have taken from the 
Woods-Saxon nuclear density. 
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To write the final formula for the cross-section we have to find the dipole distribution p(r) in 
the incoming proton, consistent with the initial condition ( 3 ). To this end we use the expression 
for the initial condition in terms of p(r) (see [3]) 

<S>(0,r,b) = ±AT A (b)g i V- 4 p(r), (5) 

where go is the strong coupling constant at a very low scale determined by the intrinsic momenta 
inside the proton. This translates into the relation in the momentum space: 

We also introduce the gluon density function 

h(y, q, b) = q 2 V 2 q( p(y, q, b) = (^-) ' <f>(y, q, b) (7) 

to finally obtain 

^(2/) = / d 2 bd 2 qf 2 (q)h(y, q, b), (8) 

oir a s0 J 

where a s o = <7o/(47r) refers to the scale inside the proton. Of course the sructure of the proton 
is unperturbative, so that a s o emerges as a new parameter in the model related to hadronic 
processes. 

Note that with these normalizations the cross-section at y = corresponds to a pure double 
gluon exchange. A better approximation consists in using an eikonalized anzatz for the initial 
condition: 

$(0,r,6) l- e -*(°' r ' 6 ). (9) 

From the formal point of view this implies including higher orders in 1/N 2 and so cannot 
be justified. Also, as follows from calculations, already at y = 1 the eikonalization becomes 
forgotten and the results with eikonalized and non-eikonalized initial conditions are practically 
the same. Still at y = they are different and the cross-section without eikonalization results 
unreasonably large. For this reason we use the eikonalized initial condition for the total hA 
cross-sections at low y and the non-eikonalized one elsewhere. 



2.2 Inclusive cross-sections and multiplicities 

In the stydy of inclusive cross-sections we have to distinguish the overall rapidity of the collision 
Y and the rapidity of the produced particle y. Our basic quantity will be the inclusive cross- 
sections lA(y,k) and lAB(y,k) to produce a jet with the transverse momentum k at rapidity 
y in a collision of a proton off a nucleus with the atomic number A or two nuclei with atomic 
numbers A and B: 

{2Tr) 2 da A{B) 

iA {B) (y,k)= dy(Pk ■ (10) 

Both cross-sections can be represented as an integral over the impact parameter b: 

lA(B)(y,k) = J d 2 bI A(B) (y,k,b). (11) 

We shall study the pA cross-section as it stands, but for the nucleus- nucleus case our study will 
be restricted to the inclusive cross-sections at a fixed impact parameter 6 = (central collisions) 
and identical nuclei, A = B. To simplify notations we shall denote the AA cross-sections also as 
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I A in all places where it cannott lead to confusion. The corresponding multiplicities at a fixed 
rapidity y will be given by 

^=hSwr Mv ' k) (12) 

for the pA case and 

» AA {y) = a/ J h = Q) J ^AAiy, k, 6 = o), (13) 

for the AA case, where a a and aAA(b) are the total inelastic cross-sections for the pA collision 
and for the collision of two identical nuclei at the fixed impact parameter b. For heavy nuclei 
one expects that <jaQ> = 0) ~ 1, so that the multiplicity (13) is just the integral of the inclusive 
cross-section over momenta. 

As argued in [4], in the perturbative QCD with a large number of colours the nucleus- 
nucleus interaction is described by a set of tree diagrams constructed with BFKL pomeron 
Green functions and triple pomeron vertices for their splitting and fusing. The structure of the 
interaction at the vertex is illustrated in Fig. 1, in which horizontal lines correspond to real 
gluons produced in the intermediate states and vertical and inclined lines describe propagating 
reggeized gluons. From this structure one sees that the produced gluons are contained in the 
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Figure 1: Interaction of three BFKL pomerons at the spliiting vertex 

intermediate states of the interacting pomerons, so that to get the inclusive cross-section one 
has to "open" these pomerons, that is to fix the momentum of one of the intermediate real 
gluons in them. For pA collisions this leaves the diagrams of the type shown in Fig 2a. For 
the nucleus-nucleus case a similar production mechanism in the old-fashioned local pomeron 
model was proven to lead to inclusive cross-sections given by a convolution of two sets of fan 
diagrams connecting the emitted particle to the two nuclei times the vertex for the emission 
(Fig. 26). The proof was based on the AGK rules appropriately adjusted for the triple pomeron 
interaction [16]. It was later shown in [17] that the AGK rules are fulfilled for interacting 
BFKL pomerons. So the same arguments as in [16] allow to demonstrate that for collisions of 
two nuclei the inclusive cross-section will be given by the same Fig. 26, that is, apart from the 
emission vertex, by the convolution of two sums of fan diagrams, constructed of BFKL pomerons 
and triple pomeron verteces, propagating from the emitted particle towards the two nuclei [6]. 
Taking into account the form of the emission vertex (see [6] ) we obtain for the pA case 

I A (y, k) = J d 2 pd 2 re ikr AV(Y -y,r)- A$(y, r, /3), (14) 



Figure 2: Typical diagrams for the inclusive cross-section in nucleus-nucleus collisions. 



and for the AA case at b = 1 

Nca. 
4 



Iaa(v, k) = I d 2 pd 2 re lkr A$(Y - y, r, f3) • A$(y, r, /?). (15) 

7T Oi r\K J 



Here ^(y, r) is a pomeron at rapidity y and of the transverse dimension r coupled to the incoming 
proton. Its normalization at y = is 

2<H0,r, b) 

(obviously the right-hand side of (16) does not depend on b). A's are two-dimensional Laplacians 
applied to ^ and $'s. 

Later from the colour dipole formalism a slightly different form for the inclusive cross-section 
was derived in [12]. For the dipole-nucleus scattering case it corresponds to changing 

2$(y, p, r) - 2$(y, /?, r) - $ 2 (y, /?, r ). (17) 

Note that in [12] it was erroneously stated that the change was from the "quark dipole" $ to 
the "gluon dipole" 2$ — <1? 2 . As seen from ( 17 ) it is not. In fact the change is equivalent to 
adding to the AGK contribution ( 15 ) a new one which has the meaning of the emission of the 
gluon from the triple pomeron vertex itself. Such a contribution is not prohibited in principle. 
From our point of view, taking into account the structure of the vertex shown in Fig. 1, its 
appearance is difficult to understand. However in this paper we do not pretend to discuss the 
validity of the two proposed formulas for the inclusive cross-sections on the fundamental level. 
Rather we shall compare the cross-sections which follow from them after numerical calculations. 

For the pA case this recipe implies taking into account a new diagram shown in Fig. 2 c. As 
a result, one finds, instead of ( 14 ), the Kovchegov-Tuchin (KT) cross-section 

IA T (V, k) = J d 2 [3d 2 re* kr A*(Y - y, r) ■ A[2d>(y, r, f3) - $ 2 (y, r, (3)], (18) 

For the nucleus-nucleus case the recipe of [12] implies taking into account two new diagrams 
for the inclusive cross-sections shown in Fig. 2 d and e. The nucleus-nucleus cross-section thus 
becomes 



1 A sligtly different coefficient as compared to [14] is due to different normalization ol $. 
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I%I(v,k) = 2 ^2 k 2 J d 2 (3d 2 re^ r [2A^(Y - y,r, f3)A4>(y,r, (3)- 



=>2/„, „ p\ A rf>2 



A$(Y - y, r, (3)A^ 2 (y, r, (3) - A<5> 2 (Y - y, r, /3)A<S>(y, r, f3)\ . (19) 

Passing in our formulas to the momentum space we find (suppressing the dependence on y 
and (3) 

A$(r) -► 2Tiq 2 A q <P{q) = 2tt%) (20) 
Introducing also a function similar to h for the pomeron 

A*(r) = 27r^°)(g) (21) 

we can express the cross-sections (14) and (15) via the gluon distributions h^°\q) and h(q) in a 
f actor ized form: 

I A (y,k) = ^^ Jd 2 pd 2 qh^(Y-y,q)h(y,k-q,f3), (22) 



and for the AA case at b = 

2 

The KT cross-section contains an additional term 



^A(y, fc) = V?T2 / d 2 (3d 2 qh(Y - y, q, (3)h(y, k - q, (3). (23) 



X(r) = -A$» = -A (27r)V> (r). (24) 
We transform it to the momentum space: 

X{q) =q 2 J d 2 re- iqr <5> 2 {r) = (2w) 2 q 2 J d 2 r{r 2 <j){r)f . (25) 



On the other hand 
So we get 

Taking into account that 



gives 



r 2 Hr) = ~J (26) 

= g 2 J d 2 q 1 A(j){q l ) ■ A4>{q - qi ). (27) 

A<f>(q) = -2irS(q) + ^ (28) 

/j2 

gl - h( qi )h(q- qi ). (29) 
9iW-^i) 2 

This should be added to the first term in the 2nd factor of the integrand in (18), which 
according to (20) is just Airh(q). So in the end in the KT cross-section instead of Airh(q) will 
appear the second term in ( 29 ). Introducing the function 

Wiq)= 2^J q^q - qi y Hqi)h{q ~ Ql) (30) 

we finally find that the KT cross-section is obtained from the AGK cross-section by the substi- 
tution h(q) — ► w(q)/2: 

I A T (v,k) = ^^ J d 2 f3d 2 qhW(Y-y,q)w(y,k-q,[3), (31) 
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Now we pass to the nucleus-nucleus cross-section. For different colliding nuclei the KT rule 
means (see Fig. 2) 

2A$ j4 (r) • A<J» B (r) -► 2A$ A (r) • A$ B (r) - A$ A (r) • A$|(r) - A$ 2 A (r) ■ A$ B (r) (32) 

where we again suppress all other arguments in $ evident from (19). Going to the momentum 
space, we have 

A$(r) -► 2?r%), -A$ 2 (r) -► -4tt%) + 2irw(q) (33) 
and we shall find the KT cross-section at a given impact parameter (5 as 

lAB(y,k,P) = 2jr2 C a 2 S k 2 J d 2 bd 2 q(h A (Y -y,q,0 -b)w B (y,k - q,b) 

+w A (Y -y,q,P- b)h B {y, k - q,b) - 2h A {Y -y,q,P~ b)h B {y, k - q, 6)) . (34) 
For central collisions of identical nuclei at mid-rapidity this simplifies to 

I%Z(Y/2, k, 0) = { d 2 bd 2 qh(Y/2, q, b) (w(Y/2, k - q,b) - h(Y/2, k - q, b)) . (35) 

So in the end the KT cross-sections can also be expressed via function h(q). 
Function h(y, k, [3) has a normalization property [6] 

j^h(y,k,P) = l (36) 

and at sufficiently high y acquires a scaling property 

h(y,kP) = h(k/Q s (y,pj), (37) 

where Q s (y,(3) is the above-mentioned saturation momentum. From ( 36 ) and ( 37 ) one easily 
establishes some properties of the new function wa- Obviously it scales with the same saturation 
momentum when h A does 

w(y,kP)=w(k/Q 8 (y,pj). (38) 

At k — > oo it has the asymptotic 

w{y,k,(3) k ^ 00 ~2h(y,k,P) (39) 

and finally 

J d 2 kw(y, k,P)=2 J d 2 kh(y, k, /?). (40) 

These properties immediately allow to make some preliminary comparison between the cross- 
sections given by ( 22 ) and ( 23 ), on the one hand, and ( 31 ) and ( 34 ) , on the other. Obviously 
if k/Q s is large both expressions give the same cross-section due to ( 40 ). In the opposite limit 
of small k/Q s , the scaling property allows to conclude that the ratio of the two cross-sections 
is a universal constant which does not depend on y, nor on A nor on (3. Our numerical results 
confirm these predictions. 
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3 Results 

Our cross-sections depend on three parameters, o"o and two coupling constants a s at the com- 
paratively large evolution scale and a s o at the very low proton scale. As mentioned, we have 
certain control on the value of o$ from the experimental data on the proton structure function at 
low x, embodied in the parametrization (3) borrowed from [15]. However we can have only very 
general ideas about possible values for the coupling constant at the two scales involved, since 
the coupling constant is not running in the model. The common wisdom is to take a s = 0.2, 
which is considered to be the value reached by the running coupling at relatively low scales. 
In our calculations we have taken this value for the evolution equation and for the production 
vertex. However with the same value for a S Q we find the total pp and pA cross-sections at low 
energies roughly 4 times greater than the experimental data. To bring them into the physically 
reasonable range we have to double the value of a S Q. So in our claculations we have taken 
a s0 = 0.4. 



3.1 The gluon densities and saturation momentum 

The gluon density in the nucleus is , up to a numerical coefficient, given by the function h(y, k, b) 
introduced in (7) [3]: 



d[xG(x,k 2 ,b)] 2N C 



d h d?k 



h(y,k,b), y = ln- 



x 



(41) 



The role of h as the gluon density is clearly illustrated in the factorized forms of the inclusive 
cross-sections ( 22 ) and ( 23 ) for pA and AA collisions. In the KT inclusive cross-section 
for hA collisions function h(y,k,b) is however substituted by (l/2)w(y,k,b), which can also be 
considered an effective gluon density taking into account the contribution from the emission 
vertex. Both densities are well defined for fixed b. As mentioned, at y > 1 they acquire a scaling 
behaviour and depend only on the ratio k/Q s (y, b) where Q s is the saturation momentum. In the 
minimum bias pA collisions and in AB collisions at a fixed impact parameter they are smeared 
out over the nucleus transverse space. In Fig. 3 we illustrate the gluon densities h and w at b = 
for A=180 and y = 2, which, with our value for a s , corresponds to y ~ 10. As one observes, 



A=180, b=0 

y=2i/(3a s ) 




k (GeV/c) 



Figure 3: Gluon densities h(k) (the left curve) and w(k) for A=180 at b = and y = 2 

both functions have the same form, with a clear maximum at a certain value of k. This value for 
h is the value for the saturation momentum Q s [18]. For w the maximum is somewhat shifted 
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towards higher momenta and the value of {l/2)w at the maximum is lower than that for h. As 
is clear from our formulas for the inclusive cross-sections, these two properties have opposite 
effect, so that in the end the inclusive cross-sections derived from AGK and introduced by KT 
result practically equal. In Fig. 4 we show the dependence of the saturation momentum Q s on b 
for A = 180 and at y = 2. Its form closely follows the nuclear density T^(6) reaching 3.8 GeV/c 



s 2 

a 

1.5 



y=2jr/(3a s ) 




b(fm) 



Figure 4: The saturation momentum Q s as a function of b for A=180 at y = 2 
at the center and dropping to 0.25 GeV/c and lower at its periphery. 



3.2 pA cross-sections 

Passing to pA cross-sections we start with the total cross-sections (Eq. (8)). They are shown 
in Fig. 5 as a function of y for A = 9, 27, 64, 108 and 180, divided by A 2 / 3 to compare with 
purely geometric cross-sections. As one observes, the cross-sections become strictly geometric 



1.8 - 




0.2 1 1 1 1 ' ' 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 

(3a s /it)Y 

Figure 5: The total pA cross-sections divided by A 2 / 3 . Curves from top to bottom on the left 
correspond to A = 9, 27, 64, 108 and 180 

at y ~ 1.2. At lower energies they grow faster than A 2 / 3 , at higher energies they grow much 
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slowlier, approximately as A - 2 . This can be explained by the role of peripheric parts of the 
nucleus, whose contribution grows with energy due to absence of non-linear damping effects and 
is relatively greater in lighter nuclei. 

Inclusive cross-sections were calculated at the overall rapidity Y = 4, which correspond to 
the natural rapidity ~ 20. Inclusive cross-sections corresponding to Eq. (14) are shown in Figs. 
6-8. Absolute values for the inclusive cross-sections are presented in Fig.6 for A = 108 at central 
rapidity (y = Y/2). At k < Q s they are oc 1/k 2 , so that the integral factor in Eq. (14) is 
practically constant. At k > Q s the cross-sections fall more rapidly. At k 3> Q s they behave 
as ~ 1/k 3 ' 1 . However this asymptotic is only reached at very high nomenta ~ 10 5 GeV/c. 
In Fig. 7 we show the inclusive cross-sections for A = 108 at different rapidities relative to 
the central rapidity. As expected, the cross-section grow towards the projectile region, since 
the rapidity of the pure pomeron factor is growing in this region. The ^4-dependence of the 




k (GeV/c) 



Figure 6: Inclusive pA cross-sections for A = 108 at Y = 4 and y = Y/2 



Y=4it/(3a s ) 
A=108 



k (GeV/c) 



Figure 7: Inclusive pA cross-sections for A = 108 at Y = 4 and different y relative to y = Y/2. 
Curves from top to bottom correspond to y/Y = 1/8, 1/4,3/8,5/8,3/4 and 7/8 
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inclusive cross-sections is different at low momenta, below Q s and high momenta, much larger 
than Q s . This is illustrated in Fig.8 where we plot (9/^4)°' 745 lA(fc)//9(A;). It is clearly seen that 




k (GeV/c) 



Figure 8: The vl-dependence of inclusive pA cross-sections at Y = 4 and y 
bottom to top on the right correspond to A = 27, 64, 108 and 180 



Y/2. Curves from 



at low momenta the cross-sections are oc A ' 745 whereas at high momenta they grow faster with 
A, approximately as A - 9 . In any case their growth with A is damped as compared with naive 
probabilistic expectations which predict them to be oc A. 

Passing to the determination of multiplicities one has to observe certain care because of the 
properties of the perturbative QCD solution in the leading approximation embodied in Eqs. ( 
14 ) and ( 15 ). As follows from these formulas inclusive cross-sections blow up at k 2 — > 
independently of rapidity y. So the corresponding total multiplicity diverges logarithmically. 
However, the physical sense has only emission of jets with high enough transverse momenta. 
Thus one has to cut the spectrum from below by some k m i n which separates the spectrum of 
jets proper from soft gluons which are not related to jets. Inevitably the multiplicity of thus 
defined jets depends on the chosen value of k m i n . We have chosen k m i n = 2 GeV/c. In Fig.9 
we show multiplicities at Y = 2 divided by A 2 ^ as a function of the rapidity y. One observes 
that on the whole the multiplicities are roughly oc A 2 ^ . They diminish with y, except at very 
small y, which is natural, since the non-linear damping of the gluon density grows as one goes 
up along the fans. 

Finally we pass to the cross-sections obtained with the KT formula ( 18 ). In Fig. 10 we 
show ratios of these cross-sections to the ones defined by the AGK rules, Eq. ( 14 ), for Y = 4 
and y = Y/4, Y/2 and (3/4)1". These ratios turn to unity at k above Q s , as discussed in the end 
of the preceding section. Below Q s the ratios are in the region 0.72-^0.87 with little dependence 
on A. Some dependence which is left can be explained by the fact that for the BFKL pomeron 
attached to the projectile the scaling regime is not valid and that for very peripheral parts of 
the nucleus it can only be reached at rapidities well above the considered ones. Due to this very 
simple relation between the two cross-sections, all conclusions about the A -dependence drawn 
for the AGK cross-section ( 14 ) remain valid also for the KT cross-section ( 18 ). 



3.3 AA collisons 



As mentioned, we only considered central collisions of identical nuclei {A = B and b = 0). The 
inclusive cross-sections obtained from the AGK rules (Eq. (15)) are shown in Figs. 11-13. Fig. 
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Figure 9: pA multiplicities at Y = 4, divided by A 2 / 3 . Curves from top to bottom on the right 
correspond to A = 9, 27, 64, 108 and 180 

11 presents absolute values for the inclusive cross-sections for A = 108 at center rapidity and 
different overall energies corresponding to Y = 2, 4 and 8. As for pA collisions, at k < Q s 
their behaviour is totally determined by the l/k 2 factor in Eq. ( 19 ), the integral factor being 
practically independent of k. At very high momenta Ia(Y, y = Y/2, k) ~ l/k p ^ with power p{y) 
diminishing with energy. From our calculations we find that p(y) ~ 3.3, 3.0 and 2.7 at y = 2, 4 
and 8 respectively. At infinite energies p seems to tend to 2 in correspondence with Q s — > oo. To 
see the dependence on the rapidity y of produced particles, in Fig. 12 we show the cross-sections 
at different y for A = 108 and the overall rapidity Y = 4. As expected the cross-sections steadily 
diminish towards the edges, their form not changing seriously. The yl-dependence of the cross- 
sections is illustrated in Fig. 13 where we show the ratio (9/A)I^/Ig(k) at mid-rapidity for the 
overall rapidity Y = 4. As in the pA case, we observe scaling at momenta below Q s , which tells 
that at such k the cross-sections are oc A with a good precision. At momenta higher than Q s they 
grow with A faster. In our calculations we found that at very high momenta the cross-section 
grow as oc A 1 ' 1 , that is much slowlier than A 4 / 3 expected from probabilistic considerations. 

Multiplicities fi{y) are obtained in AA case from the integration of Ia(v, k) over k > k m i n = 2 
GeV/c. In Fig. 14 we show them divided by A for different y at the overall rapidity Y = 4. 
We see approximate scaling. However the multiplicities in fact grow somewhat faster than A, 
which is explained by the contribution of the high momentum tail of the spectra. The change 
of the form of the y-dependence with the overall energy is illustrated in Fig. 15 where we plot 
/j,(y)//j,(Y/2) at Y = 2,4 and 8 for A = 180. One observes that the peak of the multiplicity at 
mid-rapidity becomes narrower with the growth of energy. 

All cross-sections discussed so far were obtained from the expression (15) derived from the 
AGK rules. We finlaly discuss the AA cross-sections described by the KT formula (19). As for 
the pA case we present ratios of the latter cross-sections to the AGK cross-sections in Fig. 16. 

These ratios do not practically depend on A, so we show them for A = 180 and different y 
at the total energy corresponding to Y = 4. Their behaviour is quite similar to the pA case. 
At momenta below Q s the KT cross-sections are somewhat smaller than the AGK ones. Their 
ratios grow with y from ~ 0.6 at y = 0.5 to ~ 0.8 at y = 2 (mid-rapidity). Note that at y = 0.5 
the lower part of the gluon distribution is still quite far from its asymptotic (scaling) form. So 
the values for the cross-section at this rapidity preserve a good deal of their dependence on the 
initial conditions and are not characteristic for the fully evolved dynamics of nuclear interaction. 
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At momenta much higher than Q s all ratios tend to unity, so that the KT recipe gives the same 
as the AGK one. 

4 Conclusions 

We have calculated the inclusive cross-sections and multiplicities for gluon production in proton- 
nucleus collisions and nucleus-nucleus central collisions in the perturbative QCD hard pomeron 
approach with a large number of colors. Realistic nuclear densities were employed to account 
for the peripheral parts of the nuclei, whose contribution rapidly grows with energy due to 
smallness of unitarizing non-linear effects. In contrast to the structure functions, hadronic 
processes required introduction of a new parameter in the model, the value of the strong coupling 
cosntant at very small energies corresponding to the proton structure. 

The form of the cross-sections is found to be determined by the value of the saturation 
momentum Q s , which depends on the rapidity and nuclear density. At momenta much lower 
than Q s the spectrum is proportional to 1/k 2 . Its A dependence is close to A - 7 for pA collisions 
and linear for AA collisions at b = 0. At momenta much higher than Q s the spectrum is found 
to fall approximately as l//c 2 ' 7 " 3 ' 3 with the A-dependence as ~ A ' 9 for pA and ~ A 1 ' 1 for 
AA collisions. The multiplicities are found to be proportional to A ' 7 for pA and A for AA 
collisions. Their peak at mid-rapidity for AA collisions becomes more pronounced with the 
growth of energy. 

We also compared two different forms for the inclusive cross-section, which follow from the 
AGK rules or the dipole picture. The difference between their predictions was found to be absent 
for values of momenta larger than Q s . At momenta smaller than Q s the difference reduces to a 
universal constant factor: the dipole cross-sections are just ~ 0.7-^0.8 of the AGK cross-sections. 
With the growth of energy this factor slowly grows towards unity, so that it is not excluded that 
at infinite Y the two cross-sections totally coincide at all values of momenta. All our conclusions 
about the energy, momentum and A dependence are equally valid for both forms of the inclusive 
cross-sections. 

As mentioned in the Introduction a few more phenomenological studies of the gluon pro- 
duction in nucleus-nucleus collisions were recently made in the classical approximation to the 
colour-glass condensate model [10] and in the saturation model of [11]. In both studies quantum 
evolution was neglected, so that scaling with the saturation momentum Q s was postulated rather 
than derived. The saturation momentum thus appeared as an external parameter, whose A and 
Y dependence were chosen on general grounds and whose values were fitted to the experimental 
data at RHIC. In both models the multiplicities turned out to be proportional to the number of 
participants (modulo logarithmic dependence on A). This agrees with our results. However the 
form of the inclusive distributions in momenta found in [10] is different from ours. Its behaviour 
both at small k (~ l/\/k 2 + m 2 with m = 0.0358Qs) and at large k (~ l/k 4 ) disagrees with the 
form of the spectrum we have found. For realistic nuclei the spectrum was calculated in [10] 
only up to 6-^7 GeV/c, so it is not possible to see if any change in its A-behaviour will occur 
at higher momenta. The value of the saturation momentum and the speed of its growth with 
rapidity which we have found from the QCD pomeron model with full quantum evolution are 
larger than the fitted values in both [10] and [11]. This is no wonder in view of a very high 
value of the BFKL intercept in the leading approximation which is obtained with the value for 
the strong coupling constant at present energies. From the phenomenological point of view this 
is the main drawback of the BFKL theory. To cure it one possibly has to include higher orders 
of the perturbation expansion and the running coupling constant. Although some work in this 
direction has been done for linear evolution [19], no attempts to generalize this to non-linear 
evolution in some rigour has been made yet. As it stands, the model can pretend to describe 
the data at energies considerably above the presently achieved. One may hope that future data 
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from experiments at LHC will be more suitable for the theoretical analysis in the framework of 
the model. 
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Figure 10: Ratios of the KT inclusive pA cross-sections, Eq. ( 18 ), to the ones found from the 
AGK rules, Eq. ( 14 ) at Y = 4. Curves from top to bottom refer to A = 9, 27, 64, 108 and 180 
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Figure 11: AA inclusive cross-section at mid-rapidity (y=Y/2) for A = 180 Curves from top to 
bottom correspond to the overall rapidities Y = 8, 4 and 2 
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Figure 12: AA inclusive cross-setions at different y for A = 180 and V = 4 relative to midrapidity 
y = Y/2. Curves from top to bottom correspond to y = (3/8)Y, (1/2)Y and (1/8)Y 
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Figure 13: The ratios (9/A)I A (k)/I 9 (k) for AA collisions at Y = 4 and y = Y/2. Curves from 
bottom to top correspond to A = 27, 64, 108 and 180 




Figure 14: AA multiplicities divided by A at different y for the overall rapidity Y = 4. Curves 
from bottom to top correspond to A = 9, 27, 64, 108 and 180. 
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Figure 16: Ratios of the KT inclusive AA cross-sections, Eq. ( 19 ), to the ones found from 
the AGK rules, Eq. ( 15 ) at Y = 4 for A = 180. Curves from bottom to top refer to 
y = (1/8) Y, (1/4) Y, (3/8) Y and (1/2) Y (mid-rapidity) 



